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Abstract 

QCD one-loop corrections to the full gauge invariant set of electroweak 
diagrams describing the hadronic process e + e~ — > u d s c are computed. 
Four-jet shape variables for WW events are studied at next-to-leading 
order and the effects of QCD corrections on the determination of the W- 
mass in the hadronic channel at Lep 2 and NLC is discussed. We compare 
the exact calculation with a "naive" approach to strong radiative corrections 
which has been widely used in the literature. 
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Introduction 



The measurement of VT-mass to high precision is one of the main goals of Lep 2 
and will provide a stringent test of the Standard Model (SM) [0, 0. In fact, the 
mass of the W boson in the SM is tightly constrained and an indirect determination 
of Mw can be obtained from a global fit of all electroweak data. The fit gives || 
M w = 80.338 ±0.040lg:o?lGeV where the central value corresponds to M H = 300 GeV 
and the second error reflects the change of Mw when the Higgs mass is varied between 
60 and 1000 GeV. A disagreement between the value of Mw derived from the global fit 
and the value extracted from direct measurement would represent a major failure of the 
SM. Alternatively, an improvement in the value of the P^-mass can significantly tighten 
present bounds on the Higgs mass. These studies will be continued and improved at the 
Next Linear Collider (NLC) which is expected to reduce the error on the measurement 
of W-m&ss down to about 15 MeV [|J. 

In order to extract the desired information from WW production data, theoretical 
predictions with uncertainties smaller than those which are foreseen in the experiments 
are necessary. This requires a careful study of all radiative corrections which have to 
be brought under control. In this paper we will be concerned with QCD corrections 
which are known || || to modify the shape of the W^-mass peak and the distributions 
of others kinematical variables. Since event shape variables are often used to extract 
WW production from the background it is highly desirable to have a complete next- 
to-leading order (NLO) study of these distributions for WW events. Furthermore, 
calculations of QCD corrections to 0{a 2 a 2 s ) four-jet production have recently appeared 
0. Combining these results with a complete calculation of QCD corrections to all 
0(a 4 ) four-fermion processes it would be possible to obtain NLO predictions for any 
four-jet shape variable at Lep 2 and NLC providing new means of testing perturbative 
QCD. 

W^-pair production can result in a variety of four-fermion hadronic final states. The 
simplest one is u d s c. A gauge invariant description of e + e~ — > u d s c requires in the 
unitary gauge the eleven diagrams shown in fig. 1. This amplitude is known as CC11 
in the literature. The subset of three diagrams labeled (e) and (f) in fig. 1, in which 
both intermediate IV's can go on mass-shell, is known as CC03 and is often used for 
quick estimates of WW production. Numerically the total cross sections obtained from 
CC03 and those obtained from CC11 at Lep 2 and NLC energies differ by a few per mil 
||. Other final states to which W-pairs can decay are u u d d (c c s s) and, through 
CKM mixing, uuss(uubb, ccdd, ccbb). The additional diagrams which appear 
in these latter sets contribute at most at the per-mil level to the total cross section 
and tend to produce events which resemble very little WW or even single-lV events 
since only electroweak neutral virtual vector bosons appear in them. There are also 
four-quark final states where only neutral intemediate vector boson, including gluons, 
play a role like uucc, ddss(ddbb, ssbb),uuuu(cccc) and dddd(ssss, 
bbbb). An additional important contribution to four-jet production, indistinguishable 
from four-quark events, originates from e + e~ — > q q g g. 

In this letter we present the complete calculation of QCD corrections to CC11. 
While this is only a first step in the calculation of all 0(a 4 ) four-fermion processes at 
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NLO, this reaction includes the most important source of background, namely single W 
production, providing a natural setting for a first study of the role of QCD corrections 
to gauge invariant sets of four-quark production diagrams. 

In most instances QCD corrections have been included "naively" with the sub- 
stitution rV — * rV(l + 2/3a s /n) and multiplying the hadronic branching ratio by 
(1 + oi s /iv). This prescription is exact for CC03 when fully inclusive quantities are 
computed. However, it can only be taken as an order of magnitude estimate even for 
CC03 in the presence of cuts on the jet directions and properties, as discussed in ||. It 
is well known that differential distributions can be more sensitive to higher order cor- 
rections than total cross-sections in which virtual and real contributions tend to cancel 
to a large degree. It is therefore necessary to include higher order QCD effects into the 
predictions for WW production and decay in a way which allows to impose realistic 
cuts on the structure of the observed events. The impact of QCD corrections on the 
angular distribution of the decay products of a If and their application to on-shell 
W-pair production is discussed in ref. [jj. 

Calculation 

One-loop virtual QCD corrections to e + e~ — > u d s c are obtained by dressing all 
diagrams in fig. 1 with gluon loops. Defining suitable combinations of diagrams one 
can organize all contributions in a very modular way ||. All QCD virtual corrections to 
0(a 4 ) four-fermion processes can be computed using the resulting set of loop diagrams. 




Figure 1: Tree level diagrams for e + e — > u d s c. The dashed lines are Ws. 



The real emission contribution for e + e — > u d s c can be obtained attaching a 
gluon to the quark lines of the diagrams shown in fig. 1 in all possible positions. This 
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results in fifty-two diagrams. The required matrix element has been computed using 
the formalism presented in ref. |T(| with the help of a set of routines (PHACT) |TT| 
which generate the building blocks of the helicity amplitudes semi-automatically. 

The calculation of the virtual corrections has been performed in two different ways, 
with identical results. In the first case we have used the standard Passarino-Veltman 



1 1"2|1 reduction procedure, while in the second we have used the new techniques presented 



in [y|. The matching between real and virtual corrections has been implemented as in 



using the dipole formulae of ref. [H|. All integrations have been carried out using 



the Monte Carlo routine VEGAS 15 



An important ingredient for accurate predictions of H^-pair production is the effect 
of initial state radiation (ISR). In the absence of a calculation of all O(a) corrections 
to four-fermion processes, these effects can only be included partially. In contrast with 
Lep 1 physics a gauge invariant separation of initial and final state radiation is not pos- 
sible. Only the leading logarithmic part of ISR is gauge invariant and universal. These 
contributions can be included using structure functions. Part of the non-logarithmic 
terms have been computed for CC03 and some other final states in []TB| using the current 



splitting technique |L7[ • This corresponds to splitting the electrically neutral /-channel 
neutrino flow into two oppositely flowing charges, assigning the —1 charge to the initial 
state and the +1 charge to the final state. In this way a gauge-invariant definition 
of ISR can be given. However, there are clearly cancellations between initial and final 
state radiation whose relevance is difficult to estimate in the absence of a complete cal- 
culation. Quite recently the full calculation of all O(a) corrections to e + e~ — > W + fii , fl 
has been published jn| . Non factorizable QED corrections to CC03 have been studied 
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20| . We have only included the leading logarithmic part of ISR using the ft prescrip- 
tion in the structure functions, where (3 = ln(s/m 2 ) — 1. Beamstrahlung effects have 
been ignored. We have not included Coulomb corrections to CC03, which are known 
to have a sizable effect, particularly at threshold. They could however be introduced 
with minimal effort. 

Results 

In this section we present a number of cross sections and distributions for e + e~ — ► 
u d s c. We have used a s = .123 at all energies. ISR is included in all results. The 
width of the VT-boson is kept fixed and includes 0(a s ) corrections. 

For the Lep 2 workshop [Q] the so called ADLO/TH set of cuts have been agreed 

on: 

• the energy of a jet must be greater than 3 GeV; 

• two jets are resolved if their invariant mass is larger than 5 GeV; 

• jets can be detected in the whole solid angle. 

For the NLC a slightly different set called NLC/TH has been chosen. The NLC/TH 
set of cuts differs from the ADLO/TH set in that a minimum angle of 5° is required 
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Ej > 3 GeV, y D = 0.01 


Ej > 3 GeV, y D = 0.01, 
| M m -M w |< 10 GeV 


Ej > 3 GeV, y D = 0.01, 
| M m -M w |< 10 GeV, 
smeared 


NLO 


1.1493(4) 


0.7895(5) 


0.7758(9) 


nQCD 


1.1069(3) 


1.0545(3) 


1.0479(3) 



Table I: Cross sections in pb at \fs = 175 GeV. 



between a jet and either beam (as appropriate for the larger bunch size at the NLC and 
the corresponding larger bunch disruption at crossing) and that two jets are resolved 
if their invariant mass is larger than 10 GeV. Both set of cuts will also be referred 
to as "canonical" in the following. We have preferred a different criterion for defining 
jets which is closer to the actual practice of the experimental collaborations. For mass 
reconstruction studies we have used the Durham scheme [18] , with yo — 1. x 10~ 2 at 
Lep 2 energies. At the NLC we have adopted a smaller cut y^ — 1. x 10~ 3 in order to 
have an adequate fraction of events with at least four jets. The four-momenta of the 
particles which have to be recombined have been simply summed. If any surviving jet 
had an energy smaller than 3 GeV it was merged with the jet closest in the Durham 
metric. 

Previous studies || have shown that the differences between the total cross sections 
obtained from CC11 and those obtained with CC03 are at the per mil level. Much larger 
effects have been found in observables like the average shift of the mass reconstructed 
from the decay products from the true W-mass. 

In fig. 2 we compare the NLO spectrum of the average reconstructed W-mass with 
the naive-QCD (nQCD) result at Lep 2 energies. All events with at least four observed 
jets have been retained in fig. 2. The two candidate masses are obtained forcing first 
all remaining five-jet events to four jets, merging the two partons which are closest in 
the Durham scheme, and then selecting the two pairs which minimize 

A' M = (M m - M w f + (M m - M w f . (1) 

where MmandMi^ are the two candidate reconstructed masses and Mw is the input 
W-m&ss. 

In fig. 2 the dashed line refers to the nQCD results while those of the full NLO 
calculation are given by the continous line. The corresponding cross section are given 
in table I. Fig. 2a is obtained using only the basic set of cuts described above. Since all 
experiments restrict their analysis to a region around the expected IV-mass, we have 
studied the effect of requiring that both reconstructed masses lie within 10 GeV of the 
input mass. The result is shown in fig. 2b. Finally we have tried to take into account 
some form of experimental smearing, in order to determine whether the distortion of 
the mass distribution we observe does survive in a more realistic setting. To this aim 
we have smeared the reconstructed masses entering eq. ([I]), using a gaussian distribu- 
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tion with a 2 GeV width. This procedure gives fig. 2c. It is worth mentioning that 
simulating experimental smearing at NLO is far more complicated than at tree level. 
In order to preserve the delicate cancellations between the real emission cross section 
and the subtraction terms, both contributions to four-jet quantities must be smeared, 
event by event, by the same amount. Fig. 2 shows that at NLO the mass distribution 
is shifted towards lower masses and a long tail for rather small average masses is gener- 
ated, with a corresponding reduction of the high-mass part of the histogram. This tail 
is eliminated when only reconstructed masses in the vicinity of the expected VF-mass 
are retained. Even with this additional cut, however, the NLO distribution is clearly 
different from the nQCD one. The reduction in cross section (see Table I) shows that 
a large number of soft gluons is exchanged, at the perturbative level, between decay 
products of different Ws. In our simplified treatment of experimental uncertainties 
these differences are still visible though somewhat reduced. If we try to quantify the 
mass shift using the standard quantity: 



we obtain (AM) NLO = -0.229(1) GeV and (AM) nQCD = -0.0635(4) GeV when both 
reconstructed masses are required to lie within 10 GeV of the input mass M\y = 80.23 
GeV and no smearing is applied. 

In fig. 3 we compare the NLO spectrum of the average reconstructed P^-mass with 
the naive-QCD (nQCD) result at the NLC, using the procedure already described for 
the Lep 2 case but for a smaller minimum Durham cut yo — 1. x 10~ 3 . Because of 
the larger relative momentum of the two Ws it is less likely that partons from the 
decay of one W end up close to the decay products of the other iy-boson, therefore 
the difference between the two distributions is smaller than at Lep 2 energies. 

In fig. 4 we present the distributions at y^i =175 GeV of the following four-jet 
shape variables [ETJ: 

• the Bengtsson-Zerwas angle: xbz = ^[(Pi x P2), (P3 X P4)] (fig- 4a); 

• the Korner-Schierholz-Willrodt angle : 

®ksw = l/2{Z[(pi x p 4 ), (p 2 x p 3 )] + Z[(pi x p 3 ), (p 2 x p 4 )]} (fig. 4b); 

• the angle between the two least energetic jets; au = Z[p3,p4] (fig. 4c); 

• the (modified) Nachtmann-Reiter angle: 6* NR = Z[(pi — P2), (P3 — P4)] (fig- 4d). 

The numbering i — 1 . . . , 4 of the jet momenta p^ corresponds to energy-ordered four- 
jet configurations (E\ > E 2 > E 3 > E4). We compare the exact NLO results with the 
distributions obtained in nQCD and with the results obtained at tree level from the 
standard background reactions e + e~ — > q q g g, which is the dominant contribution, 
and e + e~ — > q ± q x q 2 q 2 . 

In all subplots of fig. 4 the full NLO results are given by the continous line and the 
nQCD prediction is given by the dashed line. The q q g g and qi gi q 2 q 2 tree level 
background distributions are given by the chain-dotted and the dotted line respectively. 
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The shape variables are computed following the procedure outlined in ref. [ 2l| where the 
Durham cluster algorithm is complemented by the EO recombination scheme, namely 
if the two particles i and j are merged the pseudo-particle which takes their place 
remains massless, with four-momentum: 

Ei + Ea 

E n ew=Ei + Ej, Pnew = ; f(Pi + Pj)- (3) 

|Pi + Pil 

Using this clustering procedure all five-jet events are converted into four-jet events, 
then each event is used in the analysis if min ij= i i4 y^ > y cut with y cut = 0.008. 

In [2^| it has been shown that standard parton shower Monte Carlo programs like 



JETSET do not reproduce well the observed distribution of four-jet shape variables at 
Lep, which are instead well described by HERWIG 5. 9 A which includes four parton ma- 
trix elements. The discrepancy can be partially explained by the different distributions 
for the q q g g and the q\ q\ qi q\ final states, which can be clearly seen in fig. 4, since 
the four quark final state is included in the parton shower programs only partially. It 
is precisely in correspondence with the peaks of the q\ q\ q2 q% distributions that the 
discrepancy between data and simulations is larger. 

From fig. 4 it is apparent that at NLO four-jet shape variables distributions are 
significantly modified with respect to leading order results, which are indistinguishable 
from the nQCD distributions. It is also evident that four-jet distributions in WW 
events are markedly different from the background distributions and can be useful 
in separating the two samples. At Lep 2 energies the Korner-Schierholz-Willrodt 
angle &ksw seems to be the most effective variable for this purpose, while the angle 
between the two least energetic jets a 34 is of little use, being almost flat over the 
whole range for all samples. From our results, namely if we assume that the tree level 
background distributions closely resemble the actual behaviour of the background^, it 
appears that the differences between CC11 distributions and the dominant q q g g 
background decrease when NLO corrections to CC11 are included. 

It should be stressed that four-jet shape variables in WW events measure the 
correlations between the hadronic decays of the two W's and therefore it should be 
explicitely checked whether existing codes, NLO calculations or parton shower Monte 
Carlo programs, successfully reproduce the experimental curves. 

The shape-variable distributions at ^/s = 500 GeV are given in fig. 5. The only 
difference with respect to the Lep 2 analysis is a smaller value for the jet separation 
parameter y cut = 0.001. In particular no minimum angle between jets and either beam 
is required. A separation based on shape-variables of WW events from the background 
seems, at first sight, to be more difficult at the NLC than at Lep 2. Only for the angle 
between the two least energetic jets 034 the signal and background distribution are 
significantly different. The former peaks in the backward direction the latter is almost 
flat. On the contrary, the Bengtsson-Zerwas angle and the Korner-Schierholz-Willrodt 
angle distribution from CC11 are almost indistinguishable from those generated from 
the q q g g background. The sensitivity of the Nachtmann-Reiter angle is similar at 



2 Nagy and Trocsanyi 0, however, find large corrections, of the order of 100% for the D parameter 
and for the acoplanarity. 
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the two energies. The signal distribution peaks at small angles, particularly at higher 
energies, while the background is flatter, with a large tail which extends to 180°. The 
distributions obtained in nQCD are closer to the full NLO results than at lower energies. 

Conclusions 

We have described the complete calculation of QCD radiative corrections to the process 
e + e~ — > u d s c which are essential in order to obtain theoretical predictions for Im- 
pair production with per mil accuracy. The amplitudes we have derived are completely 
differential, and realistic cuts can be imposed on the parton level structure of the 
observed events. We have presented the distribution of the average reconstructed W- 
mass and the distribution of several four-jet shape variables at Lep 2 and NLC energies. 
The so called naive-QCD implementation of NLO corrections fails in both instances. 
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Vs = 175 GeV CC11 ISR ADLO/TH DURHAM 




x K =(M Bi +M-2M w )/2 



Fig. 2: Average mass distribution at y/s = 175 GeV. All ADLO/TH cuts are applied. 
The continuous histogram is the exact NLO result while the dashed histogram refers 
to nQCD. The corresponding cross sections can be found in table I. 
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Fig. 3: Average mass distribution at y/s = 500 GeV. All NLC/TH cuts are applied. 
The continuous histogram is the exact NLO result while the dashed histogram refers 
to nQCD. 





Fig. 4: The full NLO results (continous line) is compared with the nQCD prediction 
(dashed line) and with the tree level background distributions from q q g g (chain- 
dotted line) and q± q~i q<i q\ (dotted line). 
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Fig. 4: Four-jet shape variables at ^fs = 175 GeV, continued. 
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Fig. 5: Four-jet shape variables at ^fs = 500 GeV. The full NLO results (continous 
line) is compared with the nQCD prediction (dashed line) and with the tree level 
background distributions from q q g g (chain-dotted line) and q± q~i q2 q<i (dotted line). 




Fig. 5: Four-jet shape variables at ^fs = 500 GeV, continued. 



